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COSMIC MATTER-ANTIMATTER ANNIHILATION 
AND THE 7 -RAY BACKGROUND SPECTRUM 

F. W. Stecker 
D. L. Morgan, Jr.* 

J. Bredekamp 

Theoretical Studies Branch 
NASA Goddard Space Flight Center 
Greenbelt, Maryland 20771 

ABSTRACT 

We present here some initial results of a detailed calculation of the cosmo- 
logical y -ray spectrum from matter -antimatter annihilation in the universe. The 
similarity of the calculated spectrum with the present observations of the 7 -ray 
background spectrum above 1 MeV suggests that such observations may be evi- 
dence of the existence of antimatter on a large scale in the universe. Quantitative 
comparison of the calculations with the existing observations indicates that the 
product of interacting matter and antimatter densities at present is n p 0 n- 0 
£ 10“ 25 cm -6 . 


*NAS-NRC Postdoctoral Resident Research Associate 
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COSMIC MATTER-ANTIMATTER ANNIHILATION 


AND THE y-RAY BACKGROUND SPECTRUM 
I. INTRODUCTION 

The question of the existence of antimatter on a cosmological scale is one 
of the most basic problems of physics and cosmology. Recently it has taken on 
added interest because of the suggested possible role of antimatter in galaxy 
formation and the evolution of the universe (Alv6n, 1965; Harrison, 1967; Omn6s, 
1969) . It has long been recognized that the most promising way to search for 
evidence of antimatter on a cosmological scale is to attempt to observe y -rays 
of cosmic origin which would be produced by the decay of neutral pions arising 
from matter -antimatter annihilation. In order to determine the annihilation- 
origin of such y-rays, one must first calculate the y-ray spectrum which would 
be produced by such annihilations so that the calculated spectrum may be com- 
pared with observational data on the y- ray background spectrum for possible 
identification. 

For the purpose of discussion, we define the terms "RG" and "CR" as fol- 
lows. An RG (rest-gas) proton (or antiproton) is one possessing negligible kinetic 
energy with respect to the surrounding gas in an astronomical system. For ex- 
ample, nuclei of the galactic interstellar gas clouds from which stars are formed 
may be considered RG nuclei. A CR (cosmic-ray) proton (or antiproton) would 
be one possessing relativistic energy. We may thus envision four situations for 
discussion in which matter and antimatter might interact: 

(1) RG matter + RG antimatter (RG-RG) 

(2) RG matter + CR antimatter (RG-CR) 
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(3) CR matter + RG antimatter (CR-RG) 

(4) CR matter + CR antimatter (CR-CR) 

Although astrophysics has revealed many surprising phenomena, we would 
be hard put to imagine a system in which the number of RG nuclei did not greatly 
outnumber the number of CR nuclei. Indeed, the evolution of such a system from 
a more balanced one would violate the second law of thermodynamics. There- 
fore, if we let the number of RG nuclei in a typical system be N and the number 
of CR nuclei in the system be of the order eN, with e < < 1, then typical inter- 
action rates would have the properties 

R(RG-CR ) R (CR-RG) 

R(RG-RG) ^ R(RG-RG) ^ 6 

R(CR-CR) R(CR-CR) 

R(CR-RG) " R(RG-CR) ~ 6 

and 

R(CR-CR) 

R(RG-KG) ^ 6 

We will thus assume that the likelihood of CR-CR interactions is negligible 
compared to the likelihood of RG-CR or CR-RG interactions that will produce 
y-rays of similar characteristics, and we may limit ourselves to the discussion 
of RG-RG and RG-CR (or CR-RG) interactions. 

II. THE y -RAY SOURCE SPECTRUM 

The y -ray source spectrum from p-p interactions of the RG-CR type is 
given by 


(1) 

(2) 

( 3 ) 


2 



( 4 ) 


Oro-crCV ?) =4”np(f)Jd%I(%.r) dE ,^(E, 


V 


«£ 


^>yd **yd ^ds (^JE S ) 


For interactions of the CR-RG type, we use the corresponding expression 


Ocr-rcCV O = J d Ep I (E p , f) £ JdE.tr. (E, |E p 


L 


(5) 


=yd^yd ^ds I ) 


In equations (4) and (5), the quantity n(r) is the number of target nucleons 
in the medium per cubic centimeter as a function of position; I(E,r) is the dif- 
ferential cosmic-ray particle flux in cm -2 -s“ 1 -sr“ 1 -GeV“ 1 . The subscript p 
stands for proton, p for antiproton, s for secondary particle produced in the col- 
lision, and d for decay mode. The production function cr s (E s |E p ) represents the 
cross section for production of secondary particles of type s and energy E s in a 
collision of primary energy E p , £ yd represents the number of y-rays produced 
in the decay mode d, R yd is the branching ratio for the decay mode d (the proba- 
bility that a secondary particle s will decay via mode d) and f dt (E |E S ) is the 
normalized distribution function representing the probability that a secondary 
particle with energy E s will decay to produce a y-ray of energy E y . For the 
purposes of discussion, we will restrict ourselves to p-p interactions. 

The situation for matter-antimatter annihilations occurring near rest is 
somewhat different. In this case, we will find it more convenient to speak of 
particle densities instead of fluxes or intensities. We thus make the transformation 
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n k (r) <v> 
47T 


(6) 





n k (r) 

477 


| dv 


vf(v) = 


where v is the relative velocity between the annihilating nucleon and antinucleon, 
which we can consider to be a thermal velocity, and f (v) is a normalized distri- 
bution function representing the distribution of relative velocities between the 
interacting nucleons. 

The y-ray source spectrum from p-p interactions at or near rest is given by 
c lRG-RG( E r -r)= n p (f)n-(r) J dvf(v)v^| dE s cr s (E s |v) 

(7) 

X L ^d R r d f d,( E y l E s ) 

d 

We therefore define the emission measure B as 

B = J dln p (f)n-(r) (8) 

It follows from equation (7) that 

I RO-RG ( E y ) = 4 ^ j dlc l( ? ) “ B ( 9 ) 


where 1 is the distance along the line of sight. 


HI. CROSS SECTIONS FOR PROTON-ANTIPROTON 
ANNIHILATIONS NEAR REST 

We now state more precisely what we mean by a "rest gas." We define a 
rest gas to be a gas of particles such that no particle in the gas has an energy 
greater than 286 MeV. This rather liberal definition of rest is sufficient to 
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insure that the only secondary particles produced in RG-RG interactions that 
yield y-rays are secondary mesons produced by nucleon-antinucleon annihilation. 
Our restriction leaves out interactions of the type 

P + p-*P + p + 7T° (10) 


since the threshold for reactions of this type is 286 MeV . 


Since the threshold for nonannihilation, inelastic, p-p interactions that pro- 
duce particles other than v ° -mesons is even greater, it follows that in the RG- 
RG case, only annihilations; i.e., reactions of the form 


p + p - bosons 


L 


neutral pions 

I — *.-y- 


y -rays 


( 11 ) 


need be considered. 


Table 1 lists the experimental cross sections for free p-p annihilation as a 
function of incident antiproton kinetic energy at accelerator energies. Also 
listed are /3 cn)S and the product of the annihilation cross section (cr A ) and cms 
velocity (/3 cmg ). It can be seen from table 1 that the product of cms velocity and 
annihilation cross section is constant over a very large energy range. For non- 
relativistic energies, 

- y/ 3 os << i) (i2) 


/? being the relative velocity of the particles. We find experimentally that 

4.8 * 10 -26 cm 2 1.4 * 10“ 15 cm 3 -s -1 

A p v 


(13) 
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At all energies above 25 MeV, we find that table 1 agrees with the theoretical 
cross section 


cr 


A 



taking 


r = 0.8/ x 10“ 13 an 

p 


(14) 


(15) 


as the proton (or antiproton) radius and using the model of Koba and Takeda 
(1958) in which the nucleon acts as a black absorbing sphere of radius r p . 

At lower energies (< 10 MeV, which correspond to annihilations of matter 
and antimatter gases at a temperature <10 11 K), coulomb forces play a dominant 
role in the matter -antimatter interaction process. At energies < 100 eV, bound 
systems of protons and antiprotons can be formed as an intermediate stage be- 
fore annihilation. The situation is similar to the case of positronium formation 
(Stecker 1969). As in that case, we find that under astrophysical conditions, the 
lifetime of the p-p system against annihilation is much shorter than its lifetime 
against breakup, so that the cross section for bound-state formation becomes 
the effective annihilation cross section. Until recently the effect of coulomb 
interactions on matter -anti matter annihilation was largely neglected and theo- 
retical cr^ss sections for bound-state formation processes had not been calcu- 
lated. However, recently Morgan and Hughes (1970) have calculated the cross 
sections for the radiative capture processes 


e“ + e* -• 11 + y 

(16) 

P + P - n p + y 

(17) 
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and the atomic processes 


H + H - FI p + (Fl e or e" + e + ) (18) 

h 2 + H - ri p + (n e or e" + e + ) + (H or p + e") (19) 

H 2 + H 2 2n p + (various combinations of 2e~, 2e + ) (20) 

H 2 +H 2 -ri p + (n e or e + + e“) + (various combinations of p, p> e + , and le~) (21) 

p + H - Il p + e + (22) 

p + H 2 - Fl p + (n e or e + + e“) + (H or p + e + ) (23) 

e“ + H - ri e + p (24) 


and 


e + + H - Il e + p 


(25) 


where we have used the symbols n and Fl e to denote the bound states of the 
p-p and e“-e + systems, respectively. Reactions (24) and (25) have already been 
discussed with regard to the annihilation of secondary cosmic-ray positrons in 
a tenuous neutral hydrogen gas (Stecker 1969). The cross sections for these 
reactions were calculated by Cheshire (1964). With regard to p-p annihilation 
in intergalactic space, three processes in particular can be expected to be of 
the greatest importance: 

(1) Direct annihilations, such as (11), at energies greater than 10 MeV, 
where the cross section can be well represented by equations (14) and 
(15). 

(2) Direct annihilations at energies below 10 MeV, where the mutual cou- 
lomb attraction of the proton and antiproton distort the wave functions 
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of the two particles with a resultant increase in the direct-annihilation 
cross section. When the cross section given by equation (14) is cor- 
rected for this effect, the cross-section formula becomes modified to 


a 


A 


2TTq//3 

1 - exp (- 2 va//3) 



(26) 


where a is the fine structure constant. For p-p interactions at ener- 
gies much less than 4 MeV, where P << 2 77 a, equation (26) reduces to 


477 2 


a r 


cr. ~ 


P 2 


P « 2 77 a (27) 


Thus, at an energy of about 4 MeV, there is a transition in the velocity 
dependence of the annihilation cross section from a A ~ (T 1 to a K ~ pr 2 . 
The same is true for electron-positron interactions at energies of the 
order of 2 keV where P ~ 2 77 a. 

(3) In interactions between neutral atomic gases of hydrogen and antihydro- 
gen at thermal energies, the H-H rearrangement collision, reaction (18), 
becomes the dominant mode of p-p annihilation. Morgan and Hughes 
(1970) have calculated the annihilation cross section for this reac- 
tion and found that for energies between 10 " 3 and 1 eV (corre- 
sponding to thermal velocities at temperatures between 10 K and 
10 4 K) 

^h-IT.a - (0.31 a 0 ) 2 /T 0 M (28) 


where a 0 is the Bohr radius of the hydrogen atom. 
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Radiative capture reactions of the form of equation (17) can be shown to be 
insignificant. The cross section for annihilations of this type (Morgan and 
Hughes, 1970) is 



Equation (29) is accurate to within at least 1 percent when £ < a/7 and is exact 
in the limit a3 /3 - 0. The radiative capture cross section for the electron- 
positron system, viz, e* + e" - Fl e + y, is given by the same formula without the 
quantity (m e /m p ) 2 . For p-p annihilations at all energies of interest a A r < * cr A ; 
however, in the case of electron-po3itron annihilation (in an ionized medium) 
a A r > c; A for interactions at energies less than 20 eV. 

In figure 1 we give the croes section versus temperature, energy, and 
thermal velocity for the reactions expected to be of importance for cosmological 
problems (inter galactic p-p annihilations). Figure 2 shows the values for total- 
annihilation cross section times velocity (proportional tc the annihilation rate) 
based on t!i-^ calculations of Morgan and Hughes, which should be valid for astro- 
physical an^ cosmological problems at the temperatures and kinetic energies 
indicated.* 

*In interactions between partially ionized atomic gases of hydrogen and antihydrogen, reaction (22) 
and its charge conjugate analogue should also be considered. Morgan and Hughes find the cross 
section for this reaction to be 

- «M>57.o) J /3-' 

for temperatures between 0. 1 and 3000K. When partial ionization as a function of temperature is 
taken into account using the Saha equation, it is found that the total effective cross section is 
well approximated by (28) for temperatures £ 5000K and by (26) or (27) for greater temperatures. 
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We may now make use of equations (17), (18), (27), and (28) in order to 
reduce equation (7) to a simpler form. To do this, let us assume that mattor- 
antimatter annihilations are taking place at thermal energies in three tempera- 
ture regions defined as 

Region I: 10 11 K < T < 10 13 K 

Region II: 10 4 K < T < 10 11 K (30) 

Region HI: 10 K < T < 10 4 K 

It follows from the equations cited above that in each of these regions the 
annihilation cross section varies as some inverse power of the velocity. We can 
write an expression for this cross section for each of the temperature regions 
defined above thus 


Region I: 


cr(E s |v) ~ o-j 



Region II: 


Region ELI: 



where 

crj = 4.8 x 10“ 26 cm 2 

<7j | = 2. 2 x 10" 27 an 2 

and 

cr ni = 2.6 x iO' 18 cm 2 


(31) 


(32) 
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We will write for all three regions 


cr(E s |v) 



( 33 ) 


where Sj is the exponent of the power-law velocity dependence of a(E s I v) as 
given in equation (31) and where, as in previous discussions 

f"dE, f(E,) ■ 1 (34) 

If we assume that the velocity distribution of the interacting nucleons and 
antinucleons is thermal and can be represented by the normalized maxwellian 
distribution 


f «< v > = V^(irr) vJ exp (‘ firr) < 35 > 

at some given temperature T; and in addition, if we define a normalized energy 
spectrum of the 7 -radiation in the rest frame of the gas. 


C r f (E r0 ) * f d,< E >,: 

s J d 


(36) 


and smear it out by the distribution function, 


1/: 


D(E r lE r0 ) * 


mc‘ 

2^kT 


exp 


me 1 

2kT 


E r - E ro V 
, E r» ) 


(37) 


to take account of the thermal Doppler broadening of the 7 -ray spectrum, we 
may write equation (7) in the form 
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/ 2 V* / m \ 3/2 

^rg-rg ( E r- ? ) = n p( f ) n p(r) ^ cSi 1^“) (pfj 


J 


x I dv v ( 3 §i ) exp 


(-gf)f- 


^ ^yO ^ ( ^y l^yO ) ( 38 ) 


Equation (38) may be further reduced by elimination of the first integral to read 


q RG-RG(*V r ) n p n p a i c * ^y _ 2 




dE r0 D(E r lE r0 ) f(E r0 ) (39) 


where f (x) is the gamma function. 


In region I, 5 r = 1, and therefore, from equation (7), 

j * dv f (v) vcTj cv -1 = crj c J " dv f (v) v * v -1 = a l c dv f (v) 


- a jC 


(40) 


since J dv f (v) s 1 by definition. 

We therefore find that in region I, 

^rg-rgCV ? ) = n P ( ? ) n p(0 o-icC y | dE r0 D(E v |E y0 ) f(E r0 ) 

for 10 11 K < T < 10 13 K (41) 

In temperature regions II and HI, we cannot make the simplification given 
by equation (40); however, we can make a different simplification. In these 


12 



temperature regions, /3 << 1, and Doppler broadening due to thermal effects 
provides only a negligible distortion on the y-ray spectrum. Thus, in these 
temperature regions, we may write 


D(E r |E r0 ) - S(E y - E y0 ) 03 «1) 


( 42 ) 


and therefore 


/ dE r o D ( E r lE r o) f < E ro> - f ( E r ) («) 

In region II, 8 U = 2, and by reducing equation (39) and making use of the 
approximation (43), we obtain 

/ 2 2 V* 

%G-RG ( E r^ = n p(r) n p(r) O-II f ( E r> 

foi 10 4 K < T < 10 11 K (44) 


In region III, 8 m = 0.64, and again reducing equation (39) and using the 
approximation (43), we find 


Qrg-rg(®v r ) 


n „(r) n p(r) cr lu c£ y 


4 r (1.68) 

77 '* 


/2kT\°- 18 




me* 


f(E y ) 


for 10 K < T < 10 4 K (45) 


IV. MESON AND y-RAY PRODUCTION IN p-p ANNIHILATIONS 

Both the experimental data and the simple statistical model of Matsuda 
(1966) indicate that in p-p annihilations at rest, p -meson production is an order 
of magnitude less important than pion production and that production of other 
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mesons is at least two orders of magnitude less frequent than pion production, 
(see table 2) Table 3 shows the decay schemes of these mesons that lead to 
final-state y-rays and other relevant data. Table 4 shows some recent data 
on meson production indicating that about 20 percent of the y-rays produced 
arise through nonpionic meson production. The largest nonpion contribution to 
the y-ray spectrum is due to the p-neson decay schemes 

p* - 7T* 4- 77° ( 46) 

' — ►y + y 

The p meson is an isospin triplet (T = 1} constructed from two pions, each 
having T = 1. Evaluation of the Clebsch-Gordon coefficients for this construc- 
tion yields 

| p + ) - ( |t7 + > |tt 0 > - jrr 0 >|tt + > ) 

I/O 0 > = ( I?? 4 ) \ tt ~ > - >k + » (47) 

I/O - > = -i (k 0 > |tt- > - |t t" >|7r° >) . 

Since the p° construction does not contain any |7 t° > |7t° > terms, 

| < tt°v° |p° > | 2 = 0 (48) 

and, therefore, 

P° / TT° + 7 T° (49) 

Gamma-rays from p 1 decay (reaction (46) ) possess an average energy 
of 210 MeV, not much different from the 190-MeV average energy given to 
y -rays from the directly produced pions that are an order of magnitude more 
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frequent. Because other mesons are produced even less frequently than the p 
mesons, we can conclude that mesons other than pions have a negligible effect 
on the total y-ray spectrum from p-p annihilation. 

If we assume that the large majority of y-rays from RG-RG interactions 
arise through tt° decay, equation (36) reduces to 


f(E y ) 


I 


GO 

E>+(m£/4Ey) 




fA(S») 
(E„ 2 - "V 2 )* 


(50) 


The process 

p + p -> 77 O (51) 

is, of course, forbidden by conservation of momentum. The process 

P + P — 7T° +77° (52) 

is also forbidden since p-p annihilations at rest occur predominantly from the 
S states of the p-p system. 

The selection rule that forbids reaction (52) follows from conservation of 
G conjugation parity. 

The process of G conjugation is an extension of charge (C) conjugation, 
which holds for neutral as well as charged particles. It is defined as 

G = Ce i7rTj (53) 


where C is the charge-conjugation operator and T 2 is the second component of 
the isospin vector. From (53) we can show the commutation relation 

[T, G] = 0 (54) 
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Thus, we may describe particle states as simultaneous eigenstates of both 
G and T. 

It can be shown (Sakurai, 1964) that systems having baryon number 0 are 
in an eigenstate of G. The p-p system is just such a system. For this system 

G = (-1)L+s+t (55) 

where L, S, and T are the orbital, spin, and isospin quantum numbers of the 
state, respectively. 

In a state consisting of a single pion, L = 0,S = 0,T = 1, and, therefore, 

G = - 1. In a final state consisting of pions, G is given by 

G = (-l) u (56) 

The selection rules (55) and (56) indicate that for an S-state annihilation 
( L = 0) , the final state consisting of two neutral pions is strictly forbidden 
(Lee and Yang, 1956) . 

Therefore, the extremum y - ray energies that we would expect from this 
decay would result from interactions of the type 

p + p -• 77 + + 7T“ + 7T° (57) 

resulting in pions with the given maximum energy 

E n0, m ., = 2(2m p ) [( 2m p) 2 + m n0 ~ < 2n \,±) 2 ] ( 58 > 

3 m rr 2 

= m - — ~ 923 MeV 

p 4 m 
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Thus, the annihilation y-rays are limited to the energy region 


"2 [ E „o. 


- (E 2 - m 2 ) M 1 < E < — TE 0 

' >rO,max tt / J — y ~ 2 L ",i 


(E 2 

nO, max 


-"V 2 )*] (59) 


or 5 MeV < Y. y < 919 MeV . 


The process 


p + p -* y + y 


(60) 


may also occur; but this process, involving the electromagnetic emission (em) 
of two photons, is of the order 


2.6 x 


10 -30 cm 2 


- a cr 


A, s t rong ” 


(61) 


and therefore this process may be neglected as a significant contribution to 
y -ray production. 


V. THE GAMMA -RAY SPECTRUM FROM 

PROTON-ANTIPROTON ANNIHILATIONS AT REST 

As we showed in the previous section, we can neglect the contribution from 
the decay of mesons other than neutral pions in calculations of the y-ray spec- 
trum from p-p annihilations at rest. The normalized y-ray spectrum was cal- 
culated numerically from the relation (50) with the normalized distribution 
function f A (E n ) taken from the calculations of Maksimenko (1958), based on the 
statistical theory of multiple particle production. 

The resultant spectrum, up to 750 MeV, is shown in figure 3. Frye and 
Smith (1966) have recently calculated the y-ray spectrum from p-p annihilation 
up to 500 MeV, based on recent measurements of charged pions from p-p anni- 
hilation by the Columbia University group. The excellent agreement between 
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the results of figure 3 and the calculations of Frye and Smith not only serves as 
a mutual check on the calculations but also supports our previous conclusion 
that mesons other than pions have a negligible effect on the total y- ray spectrum 
from p-p annihilation at rest. 

VI. COSMOLOGICAL REDSHIFTED SPECTRA 
Our calculations are in accord with the suggestions proposed by the models 
of Harrison 2 and Omnes 3 . These models assume that tne universe consists of 
equal amounts of matter and antimatter which were separated info distinct re- 
gions at the earliest stages in the big-bang model of the universe. In the appen- 
dix it is shown that in the energy region of observational interest discussed 
here (E y > 1 MeV), where the redshift of origin of the y- ray background is < 100, 
we can assume that annihilations take place on the boundaries of colliding regions 
of matter and antimatter so that the total annihilation rate is proportional tc 
(1 + z) 6 . We use the commonly defined cosmological parameters H 0 and f 1 
where Hg is the Hubble constant and Q is the ratio of the average atomic matter 
density in the universe to the critical density n c ~ 10 * 5 cm -3 needed to gravita- 
tionally close the universe. We also define the quantity £ which denotes the 
mean fraction of the total atomic density interacting at the boundaries between 

the regions of matter and antimatter. The atomic densities n and n_ are each 

p p 

proportional to (1 + z) 3 and it is assumed that on the average n p = n_ = Qn,.. 
Thus we can define a mean density of interacting matter (or antimatter) n int 
given by n ,nt = £fin c . 

For cosmological y - ray production, the quantity B defined by equation (8) 
then becomes 

, f . dl 

B - I di n p (z) n-(z) (62) 
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z is defined as the redshift suffered by a photon originating at a distance l. 


The derivative 


dz H 0 (l + z) 2 (1 + fiz)* 


n P ( 2 ) = Hp.oC 1 + z ) 3 


n p( z ) = n p, 0 ( l + z )“ 


(see, e.g., McVittie 1965). 


Defining 


B 0 * £ 2 


( H „) n p.o n P. 


we find by taking account of the effects of redshift, expansion, and time dilatation 
in the expanding universe 


b o ct . e Sl /2k\(* »') / ®i \ 

Ia(E > ) 2 tt 3/2 \ m ) ( 2 2/ 


X dz [T(z)] 


( 2 ) ±i * z _ ) _ G A [(1 + z) E l (65) 

(1 + fiz)* AL J yJ ' 


where i = II, HI; cr n = 2; 8 ln = 0.64; cr ll = 1.1 x 10~ 27 cm 2 ; cr m = 2.6 x 10“ 18 
cm 2 ; G a (E y ) s f k (E y ) where f A (E y ) is shown in figure 3; and T m is the matter 
temperature of the universe, which is a function of z. The matter temperature 




f 2. 7K (1 + z) for z > 150 


T( z ) 


L 1. 8 x 10- 2 K (1 + z ) 2 for z < 150 


( 66 ) 


(Zel'dovich, et al., 1969). 

It follows that equations (65) and (66) together with the bounded form of the 
function G A (E y ) yield source spectra of the form E' 1 as shown in table 5 when 
absorption effects can be neglected. 

At high redshifts, when pair production and Compton scattering become 
important, it becomes necessary to solve a cosmological -photon-transport 
(CPT) equation in order to determine the y-ray spectrum. 


VII. THE COSMOLOGICAL-PHOTON-TRANSPORT EQUATION 

The cosmological-photon-transport (CPT) equation for x-rays and y-rays 
has been discussed at great length by Arons (1971a, b). Our formulation differs 
from his basically only in the definition of the y-ray intensity, I. Our intensity 
is defined as a photon intensity in units of photons per cm 2 -s-sr-MeV (or in 
energy units of 0.511 MeV = m e c 2 ). The intensity defined by Arons is in terms 
of an energy flux rather than a photon flux. The energy flux decreases more 
sharply with redshift than the photon flux, since the energy of each photon suffers 
an adiabatic expansion loss as the universe expands. Thus, the form of the CPT 
equation used here is slightly different than that used by Arons. 

Defining y ■ 1 + z, and e = E y /m e c 2 , we can write the CPT equation as 
follows 
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where y -3 I is the specific comoving photon intensity in the expanding universe, 
Q A (e,y) is the source term of y -rays produced in matter -anti matter annihila- 
tions, and n(y) is the average gas density in the intergalactic medium. The 
quantity de/dt describes the energy loss rate of photons due to the expansion 
of the universe which is given by 


de 

dt 


e dy 
y dt 



H 0 y 2 ^ 1 + 0 (y - 1) 


( 68 ) 


(see, e.g., McVittie 1965). The second term on the right hand side of equation (67) 
describes the absorption of > -rays by intergalactic gas due to pair -production 
and Compton scattering. The integral term describes the transport of photons 
being scattered down in energy by Compton interactions. 

The partial derivative 

3 dy 3 « i 3 

SV'd!^ = - H„y J V 1 + n(y- 1) (69) 


Substitution of equations (68) and (69) in equation (67) and subsequent simpli- 
fication yields the CPT equation in the form 
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where 


H 0 

cr = ~ 10" 2 3 cm 2 

u n c 

and 


(71) 


n c = 3H 0 2 /8 ttG ~ 10" 5 cm" 3 . 


G being the Gravitational constant. 


The upper limit on the Compton-scattering integral, 


b(e) 


e/(l - 2 e), e < -z 


00 ' e ~ 2 


(72) 


and the interaction velocity function 
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and G (e) is the 7 -ray source term normalized per interaction as shown in 
figure (3) . 


We further define 


r- 


V = 


= 2.5 x 10" 2 


(74) 


where r Q is the classical radius of the electron. We then introduce a dimension- 
less absorption term A (e) = A c (e )+ A p (e) to take account of both Compton 
scattering and pair production where 
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Vffl. COMPARISON OF CALCULATED 7 -RAY SPECTRUM FROM 
COSMOLOGICAL MATTER -ANTIMATTER ANNIHILATION 
WITH C '"SERVATIONS 

The function I A (E y , y = l)/n p L n , obtained by numerical solution of the 
CPT equation using the method of characteristics, is shown by the solid line in 
figure 4. If we chose the values fl = I and £ = 3 x 10" 8 , we obtain a theoretical 
annihilation spectrum compatible in both form and intensity to the observed 
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spectrum of cosmic background 7 -radiation above 1 MeV as is shown in curve A 
of figure 5. We conclude that tne observations of Vette, et al. (1970), which indi- 
cate a marked deviation from the power-law X-ray spectrum (curve X), may be 
evidence of the existence of antimatter on a cosmological scale. This conclusion, 
however, is subject to the following conditions: 

(1) There must be an additional power-law background component of X-rays 
below 1 MeV of another origin. Our calculations of the annihilation 7 -ray spec- 
trum below 1 MeV (not presented in their entirety here) fall well below the ob- 
served power-law spectrum in that energy range 

(2) The peak in the calculated annihilation spectrum near 1 MeV is caused 
by absorption and scattering of the 7 -rays by interactions with an intergalactic 
medium having an average density near the critical value n c ^ 10“ 5 cm' 3 . 

(3) There is strong evidence that at a redshift of z ~ 2.5 (or perhaps some- 
what greater) the intergalactic medium is strongly ionized and equation (3) does 
not hold (Gunn and Peterson 1965, Rees 1969). At the redshift where reionization 
occurs there will be a rapid drop in the cross section for annihilation (see 
figure 1 ) resulting in a sharp steepening of the 7 -ray spectrum at an energy 

E y c ~ 50 MeV. Above E y c the spectrum will drop off much more rapidly than 
the power-law dependence ~E “ 2 - 86 shown in figure 3. Assuming that the re- 
ionization of the intergalactic medium takes place at a redshift of 2.5, we obtain 
the dashed line shown in figure 5 for the high-energy end of the annihilation 
spectrum. 

(4) The observational data on the 7 -ray spectrum above 1 MeV do not 
provide enough information at present to yield a unique identification of their 


24 


% 

B- ' 

fV 

I 

■g 





ft; 


v 

§./• 


S. 



- 


origin. There have been several alternative attempts to explain the background 
flux as being due to the decay of neutral pions produced by cosmic-ray inter- 
actions at high redshifts (the protar hypothesis (Stecker 1969a, 1969b, 1971)), 
galactic electron bremsstrahlung (Rees and Silk 1970), extragalactic electron 
bremsstrahlung (Silk 1970) extragalactic proton bremsstrahlung (Brown 1970) 
and nuclear emission lines (Clayton and Silk 1969). More recent work shows 
that none of these alternatives except the protar hypothesis seems capable of 
explaining both the form and intensity of the observed flux (Vette, et al. 1970, 
Stecker, et al. 1971, Stecker and Morgan 1971, Jones 1971).* The form of the 

-ray spectrum produced as a result of cosmic-ray interactions at high redshifts 
is shown in the curve marked CR in figure 5. It can be seen from figure 5 that 
definitive measurements of the cosmic 7 -ray background spectrum above 6 MeV 
will provide an observational test between the protar hypothesis (curve CR) and 
the annihilation hypothesis (curve A). The annihilation hypothesis predicts a 
steepening in the background spectrum above 50 MeV whereas the protar hy- 
pothesis predicts a steepening in the spectrum at about 7 GeV (Fazio and 
Stecker 1970). Should future investigations indicate that the background 
7 -radiation above 1 MeV is not due primarily to matter-antimatter annihilation, 
this will place an upper limit on the quantity £ of 3 x 10 -8 or alternatively an 
upper limit on the product of interacting matter and antimatter densities at 
present of n p in o n^"^ ~ 10 -25 cm -6 . 


*A suggested explanation by Sunyaev (ZhETF Pis. Red. 12, 381 (1970)) has recently come to our 
attention. He suggests that bremsstrahlung radiation from relativistic electrons in the nuclei of 
Seyfert galaxies would have a spectrum consistent with the Vette, et al. results and have a sharp 
cutoff above 20 MeV consistent with the upper limits at higher energies. The sharp cutoff at 20 
MeV suggested by Sunyaev has been disputed by Prilutsky, Ochelov, Rozental and Shukalov 
(U.S.S.R. Academy of Sciences Institute for Space Research Preprint #51, 1970). The problem of 
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producing a high enough intensity of gamma-rays from this mechanism has been discussed by 
Bisnovaty-Kogan and Sunyaev (Ap. Lett. 7, 237 (1971)). 

Note on recent experimental results: Trombka (private communication) has indicated that a recent 
reanalysis of the Ranger 3 results above 1 MeV is in agreement with the measurements of Vette, 
et al. Golenetsky (Moscow Seminar on Cosmic-Rays and Astrophysics 1971) has reported that new 
upper limits obtained by a detector aboard the low-altitude Cosmos 135 satellite for "y-rays of 
energies up to 3 MeV are in conflict with the measurements of Vette, et al. Daniel (l2th Int'l. 

Conf. on Cosmic Rays, Hobart, Australia ^August, 1971, OG-27) has reported 
fen measurements between 0.25 and 4.2 MeV which appear to be consistent with 
the measurements of Vette, et al. Measurements made aboard Apollo lp are 
now being analysed and may be available in the near future. Uncertainties 
in the present data indicate the need for cleaner, more sophisticated 
measurements in the future. 
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Table 1 


Experimental Cross Sections for p-p Annihilation as a 
Function of Incident Antiproton Kinetic Energy 


T MeV 

a K' 

mb 

/j_, cms 

/3_,ct a , mb 

Reference 

25 to 40 

192 

± 34 

0.13 

25 

Loken and Derrick (1963) . 

45 

175 

± 45 

.15 

26 

Cork et al. ( 1962) . 

40 to 55 

155 

± 27 

.16 

25 

Loken and Derrick ( 1963) . 

55 to 80 

118 

± 26 

.20 

24 

Loken and Derrick (1963) . 

90 

101 

± 9 

.22 

22 

Cork et al. (1962) . 

145 

99 

t 8 

.28 

28 

Cork et al. ( 1962) . 

245 

66 

± 6 

.36 

24 

Cork et al. (1962) . 

7000 

23.6 ± 3.4 

~ 1 

24 

Ferbel et al. (1965) . 


Table 2 


Average Energies and Production Rates for Various Particles 
Produced in p-p Annihilations at Rest 
[Matsuda, 1966] 


T, 

particle or pair 

Average 
energy, MeV 

Calculated 
production rate 

Experimental 
production rate 

TT 

380 

3.96 

3.94 ± 0.33 

P 

850 

2.3 x 10" 1 

(2.5 ± 0.6) x 10“ 1 

U i 

940 

4.2 x 10“ 2 

(4.5 ± 0.7) x 10" 2 

V 

860 

2.4 x 10“ 2 

(1.4 ± 0.5) x 10~ 2 

KK 

(K)660 

3.1 x 10“ 2 

(3.3 ± 1.6) x 10“ 2 

KK*, KK* 


20.6 x 10" 3 

(8.8 ± 1.8) x 1C -3 

K* K* 

(K*)980 

3.8 x 10‘ 3 

(3.9 ± 0.7) x 10“ 3 


31 




Table 3 


Decay Modes. Branching Ratios, and ?-Ray Multiplicities 
for Various Particles Produced in p-p Annihilations 


Decay mode 

Branching ratio, 
R 


p 4 -• 77 4 + 7 7° 

~1.00 

2.0 

CO - 77* + 77“ + 77° 

.89 

1.78 

-• 77° + y 

.10 

.30 

rj - y + y 

.386 

.78 

— 3tt° or 77° + 2 y 

.308 

1.5 

77 + + 77~ H 77® 

.250 

.5 

-• 77 + + tt" + y 

.055 

.05 

K 4 - 77 4 + 77° * 

.215 

.43 

r K0 - 77° +77° 

.155 

.62 

) K° - 77° + 77° + 77° 

.133 

.80 

K° - 77* + 77" +77° 

.067 

.1 

K* - K + 77 

— 

— 


Other y~ ray-producing decay modes have negligible branching ratios. 




Table 4 


Production Rates for Various Meson-Producing Channels 
in p-p Annihilations at Rest 
[From Baltay et al., 1964] 


Channel 

Rate, % 



O 

o 

1.4 

± 0.2 

2 

2.8 

P* Tl 4 

2.9 

x 0.4 

2 

5.8 

0 
♦ 

1 

m 

5.8 

+ 0.3 
- 1.3 

0 

0 

P° P° 

0.4 

± 0.3 

0 

0 

rr* n“ tt° 

7.3 

± 1.7 

2 

14.6 

P * T* TT* 7T~ 

6.4 

± 1.8 

2 

12.8 

J> TT* 7T~ * 

3.8 

± 0.4 

2 

7.6 

7J° V* TT~ “ 

1.2 

± 0.3 

~ 2.6 

3.1 

00° p° 

0.7 

± 0.3 

A 

£* 

— 

© 

o 

0.22 ± 0.17 

•'•2.6 

— 

“ Includes cases where n* v~ were 

from 

p° decay. 





Table 5 



Exponents of Cosmological Power-Law Annihilation ?-Ray Spectra 

Q z 


70 MeV 

E y0 - 

(1 + z) 

5 

r 

1 0-150 


500 keV-70 MeV 

0.64 

2.86 

150-1000 


70 keV-500 keV 

0.64 

2.68 

10 3 -10 6 


< 70 keV 

2 

2.00 

~0 0-150 


500 keV-70 MeV 

0 64 

3.36 

150-1000 


70 keV -500 keV 

0.C4 

3 18 

10 3 -10 6 


< 70 keV 

2 

2.50 
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FIGURE CAPTIONS 


Figure 1. Proton-antiproton annihilation cross sections as a function of tem- 
perature and kinetic energy for free nucleons and atoms (Morgan and 
Hughes, 1970). The curve marked a p takes into account the mutual coulomb 
attraction between the proton and antiproton; the curve marked a is an 
extrapolation from the accelerator data that fails to take this effect into 
account. 

Figure 2. Cross section times velocity for hydrogen-antihydrogen annihilation 
of nucleons given as a sum for the interactions shown in figure 1 integrate 1 
over a Maxwell-Boltzmann distribution at temperature T and taking ioniza- 
tion into account. (a) °h + h (extrapolated), (b) cr H+ „ (calculated), (c) ^ h+ h, 
o'? + — , cr_ +H and cr p+ _, multiplied by appropriate factors to take account of 
fractional ionization as determined using the Saha equation, (d) a p+ - , taking 
coulomb attraction into account, (e) cr p+ -, from high-energy accelerator 
data. 

Figure 3. Normalized y-ray spectrum from p-p annihilation. 

Figure 4. The cosmological y-ray annihilation spectrum calculated by numerical 
solution of the CPT equation (equation (80)) for fi = 1. The solid line repre- 
sents the complete solution. The other curves show the effect of neglecting 
the absorption and scattering terms in the CPT equation. 

Figure 5. Observations of the y-ray background spectrum together with the 
theoretical spectra discussed in the text. 
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TEMPERATURE (K) 


KENETIC ENERGY IN CENTER OF MASS SYSTEM (eV) 


Proton-antiproton annihilation cross sections as a function of temperature and kinetic 
energy for free nucleons and atoms (Morgan and Hughes, ’970). The curve marked a 
takes into account the mutual coulomb attraction between the proton and antiproton; P 
the cjrve marked cr , is an extrapolation from the accelerator data that fails to take 
this effect into account. 
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(PHOTONS/CM SEC Sr MeV) 


SOLUTION WITH 
TRANSPORT AND 
ABSORPTION 


\no ABSORPTION 
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\ 
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/ WITHOUT TRANSPORT 
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E y (MeV) 


Figure 4. The cosmological y*ray annihilation spectrum calculated by numerical solution ot the 
CPT equation (equation (80)) for Q = 1. The solid line represents the complete so- 
lution. The other curves show the effect of neglecting the absorption and scattering 
terms in the CPT equation. 
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APPENDIX: THE BARYON SYMMETRIC COSMOLOGY OF OMNES 


AND THE VALUE OF £ 

J. L. Puget 
F. W. Stecker 


1. The Omnes Model 

According to the calculations of Omnes, (1971a,b), in black body radiation 
at high temperature (>350 Mev) there is a phase transition which leads to the 
separation of matter and antimatter. In the big bang model of the universe this 
corresponds to a time 10“ 5 sec after the initial singularity. 

During the period when the temperature remains higher than a fraction of 
an Mev (until the end of the leptonic era) the nucleons and antinucleons diffuse. 
The annihilation is only limited by the speed of diffusion and the statistical 
fluctuations which depend upon the size of the initial separation. The next period 
is called by Omnes the coalescence period. The study of the behavior of an 
emulsion of matter and antimatter by Omnes (197 la, b) shows that the charac- 
teristic size of a region of pure matter (or antimatter) in the emulsion is inde- 
pendent of the initial mechanism and only depends on the ratio of the density of 
nucleons to the density of black body photons as a parameter. This ratio tj = 

n P / n ph is 10 ' 8 for n = !• 

The coalescence period ends at a temperature ~jeV when the plasma com- 
bines, and the size of the emulsion regions is fixed by the time at which the 
plasma forms neutral atoms and is ^k y wherr, \ y is the mean free path of the 
y rays at that time. 
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2. The Neutral Period 


Omnes (1971a) considers that the start of the neutral period in his cos- 

the 

mology occurs at the same time as impure matter big bang model. We will now 
show that this is very inaccurate because the annihilation process results in 
ionizations which delay the start of the neutral period. Peebles (1968) has given 
a detailed solution of the combination problem when the plasma is only composed 
of matter. The equation he gives for the electron density is: 

Bi-b 2 \ 

kT J C ( T ) 

where B n is the binding energy for hydrogen with the electron in an orbit corre- 
sponding to the n th principal quantum number and n H is the density of hydrogen 
atoms in the ground state. The a c ,/3 c and C (T) coefficients are functions of 
the temperature T . 

The annihilation rate per unit volume is p, = — where R is the character- 

I\ 

istic dimension and is the annihilation rate per unit surface area of the an- 
nihilation layer, 

v 

with h given by 

h 

and D is the diffusion coefficient. 

This diffusion coefficient has been discussed by Omnes (1971a). He finds 

v % - N 0.48 v t 


= h 


n P °p <<T v> 


= [ 24 idh]* 


dt n « 
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where v is the thermal velocity and N is the density of nucleons (antinucleons) 

at the boundary of the annihilation layer. N is lower than n p (or n p ) because of 

the effect of th. e annihilation pressure (Omnes 1971b). During the coalescence 
N 

period the ratio — is almost constant and is ~ 0.05. 
n 

p 

In order to investigate the effect of this annihilation on the recombination of 
the plasma one must consider the behavior of all the products of the annihilation. 
The neutrinos escape carrying ~ 50% of the energy. They -rays lose their 
energy by Compton collisions with the electrons of the medium with an effective 
cross section <cr y > ^ ^ cr T (a T is the Thompson cross section). The secondary 
electrons have an energy of the same order of magnitude as the y-rays. y-rays 
are also produced by annihilation of the electrons and positrons after they have 
been slowed down. 

The electrons and positrons undergo two different types of interactions; 
inverse Compton interactions on the black body radiation with a cross section 
cr T , and collisions with other electrons (free or bound). For a value of the pa- 
rameter t) = 10" 8 and for z ~ 1000 the first process is dominant for relativistic 
electrons. Using a simple representation of the spectrum of the electrons pro- 
duced in the annihilation; 0.86 electron with an energy 40 MeV, 0.41 with an 
energy 135 MeV and 0.1 with an energy 300 MeV, one gets a rough idea of the 
number and the energy of the X-rays produced for one annihilation: 10 4 with 
the energy 20 keV, 2.10 4 with the energy 5 keV and 5.10 4 with the energy 0.5 
keV . This estimate is based on the assumption that the average energy of the 
X-rays produced in a black body radiation at the temperature T by relativistic 
electrons of energy E is 
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(e.g. Ginzburg 1969) 


e 


X 


k T 



2 


These X-rays can themselves undergo two interaction processes, photo- 
ionization and Compton collisions with free electrons. The photoionization cross 
section depends strongly upon e x 



7/2 


(see, e.g., Heitler 1954) 

These X-rays are very efficient ionizing agents and tend to slow down the 
recombination of the plasma. At the beginning of their life, the most energetic 
X-rays are degraded by Compton collisions with the free electrons. This 
process becomes unimportant when the energy of the X-rays falls below a few 
keV at which point photoionization becomes important if only one percent of the 
hydrogen is neutral. In this way, half of the energy released in the annihilation, 
which is carried by the ?-rays and the electrons goes into thermal electrons in 
a short time compared to the cosmological time. 

The number of ionizations per annihilation is thus given by 

m c 2 (eV) n H 

k i = ^ 

1 32 n 

p 

because one ionization corresponds to a loss of 32 eV by a fast electron. 

The equation for the recombination equilibrium then becomes 

^ n e = a c n,2 C(T) - p m • kj = 0 
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which reduces to 


C(T, x) x 2 = 2.5 10' 2 — (1 - x) 

10 4 

where 

n e 

X = — . 

n 

It follows from Table 1 that the recombination is very slow compared to the 
big-bang model without matter-antimatter symmetry. It follows that the final 
size of the matter and antimatter regions is correspondingly larger than that 
calculated by Omnes (1971a). 

However, two other effects must be considered in determining the final size 
of these regions: 

1. E. Schatzman (1970) has shown that even very small magnetic fields 

can drastically alter the annihilation rate. When matter and antimatter annihilate 
on an interface a natural mechanism leads to the building of a magnetic field 
(see below). 

2. If the plasma remains ionized much longer than Omnes (1971a) supposed, 
it is not certain that the annihilation provides enough momentum for the coa- 
lescence to continue until neutralization is complete. 

We will not deal here with these problems which we hope to discuss later. 
Omnes (19 7 la) concluded that the final mass of these regions would be of the 
order of large galaxies. We conclude that it must be larger because of the 
above discussion and can be of the order of galaxy clusters. 
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This conclusion is supported by the fact that the collisions of galaxies and 
antigalaxies in dense clusters, if we suppose their number equal, leads to a 
y -ray flux above 100 MeV larger than the Clark et al. (1968) limit by a factor 
10 3 , (Puget J971, unpublished). 

3. Metagalactic Annihilation As a Function of the Red-Shift 

3.1 Initial Conditions 

At z somewhere between 300 and 1000 the effect of neutralization and of 
annihilation pressure is to stop the coalescence by separating matter and anti- 
matter into independent clouds. We suppose that these clouds are the proto- 
clusters of galaxies. Wc know their actual number, 7l 0 , is ~ 5 x 10' 5 Mpc" 3 
(Abell 1965). Denoting the red-shift corresponding to the time of separation of 
the clouds as a parameter z i , we have computed the radius of the clouds at this 
time, and their corresponding mass as determined by the relation 

3 ff R 3 \ ( 1 + z i ) 3 = 1 

where R(Zj)~ (z t ) is the final radial dimension of the separate regions as 
previously discussed. The results are shown in Table 2. 

3.2 Cloud- Anti cloud Collisions 

For z < 200 the cloudf e well separated and the annihilation rate is di- 
rectly related to the collision between clouds and anticlouds. When a cloud of 
neutral matter collides with a cloud of neutral antimatter, all the atoms from 
one cloud reaching the other cloud are annihilated; the rate of annihilation is 
given by = n H v (v is the relative speed). The strong flux of electrons and 
>-rays produces by Compton interactions a current of the same sign on both 
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sides of the annihilation layer and one can see that, for the collision of two 
roughly spherical clouds, this current produces a very low primary magnetic 
field but one which is strong enough so that the Larmor frequency is greater 
than the inverse electron lifetime by Compton losses. (Puget 1971, unpublished) 
Under these conditions, the current is not normal to the annihilation layer but 
parallel to it and we have to consider the behaviour of the magnetic field in re- 
lation to this drift current as the ionization front propagates. The idea of the 
following comutations is due to E. Schatzman (1970). The general process 
remains the same but some details have been corrected. 

Figure 1 depicts the processes involved. Two ionization fronts are moving 
symmetrically with speed c on both sides of the annihilation layer. 

This ionization is mainly due to photoionization by the X-rays. The elec- 
trons ejected in this process have an initial energy '~e x and produce a drift 
current when they are slowed down. All the characteristic times for these 
processes are shorter than r, the lifetime of the primary electrons produced 
in the annihilation. ( Tis also the characteristic time for the production of the 
X-rays.) 

Between the annihilation layer and the ionization front, a flux of X-rays 
and a flux of y-rays produ. last electrons by Compton collisions. These fast 
electrons, when slowed down, generate a drift current J of opposite sign on each 
side of the annihilation layer. 

per unit volume 

In the plasma, the particles are in equilibrium between the force/F x due to 
the X- and 7 -rays and the electromagnetic force ptr nit volume 
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Maxwell's equations together with the fact that all the field quantities are 
independent of y due to the choice of the frame of reference (indicated in Fig. 1) 
yield the relation 


BBz v x 3Bz 

TT + = 4 ” F 

Z 


(A-l) 


if we suppose that the drift velocity of the plasma, v x , is almost independent 
of t. 

The force F x has two components coming from the X-rays ana the y-rays. 


n e n e 

= L “.ST 'VS + E v . CT rT n . 


where n y is the number of y-rays of energy e y produced in one annihilation. 
And where we have used a simple representation of the annihilation y-ray 
spectrum assuming 

50% production with an energy of 70 MeV 

35% production with an energy of 250 MeV 
and 15% production with an energy of 500 MeV , 
to estimate the value of F x . 

The solution of the Equation (A-l) gives 

B 2 

,ar cF x (x + v x t) for x < ct - cr 


and 


& 

877 


& 



(ct-c-r) 


+ f;v x t [1 - * 2 ] 


for ct 


CT < X < Ct 
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where 


- x + c t 


and is the force in the ionization front, mainly due to photoionization pro- 
cess of the hydrogen by the X-rays, 

r 2e * 

r = J <t> n„ cr, — d £ x 

where 4 > is the differential flux of X-rays. 

B» 

The quantity — as a function of x is given on Figure 2. 

The time necessary for developing a magnetic field strong enough so that 

B 2 > 1 2 

Frr - 2 n V 

is given by 

t _ 1 2 1 

t B - -j n m p v 2 

This time is muc* larger than r and justifies the approximations made. 
The field and the charged particles move together. When the energy density in 
the field in one cloud is much larger than the kinetic energy of the other cloud, 
the field cannot propagate to the other side of the annihilation layer. The an 
nihilation stops and the collision becomes almost elastic with a ’’magnetic 
cushion" between the two clouds. The surface area of the two clouds in contact 
is a function of the impact parameter, b, and the radius, R given by 
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s ~ 


(A-2) 


tt( 2R - b V b 2 \ 

—2— [ R2 -t) 

where R is the radius of the clouds. The annihilated mass is 

1 0 0 , 1 
M a = 2 n m P v S TV 


The fraction 


i s 



1 m 2 v 2 
u P 

M t 8X10" 28 


(A— 3) 


The number of these clouds per unit volume is ft = ft Q (1 + z) 3 . ft Q is their actual 
density can be evaluated as: 5 10“ 5 > ft > 10“ 5 (Mpc“ 3 ) depending on the value 
of the Hubble constant. 


The average density of matter in the universe is given by 


n( z ) - n 0 (! + z ) 3 

so the mass of the clouds is fixed by the value of ft 0 and the hypothesis 0 = 1. 
The total mass M t = 2.8 10 48 g and 1.4 10 49 g with ft 0 given by 5 x 10* s and 
10" 5 Mpc“ 3 respectively. The relative cloud velocity v is chosen to be 3 x 10 8 
cm/s (Abell, 1961). During most of the collision period, R is almost constant 
because of the strong cooling by the black body radiation. R is certainly in- 
creasing in a non-negligible manner for low z (< 10) but other changes appear 
at the same time: strong decrease of H - H cross-section (Bardsley 1971), 
formation of galaxies, existence of a magnetic field created in the previous 
collisions and very low probability for intercloud collision. We conclude that 
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the behaviour of the annihilation rate at low red-shift is not clear and the 
formula given below can be wrong by a large factor at very low red-shift. 

3.3 Annihilation Rate as a Function of z. 

The number of collisions is given by: 

dN - ^ Jl 2 27rbdb v 

where b is the impact parameter. 

The annihilation rate is: 

f 2R ft 2 M a 

P. = 1 2nbdb v - (A-4) 


where M a = M a [*(b) ] . 

Using Equation (A-3) , we can define 


<i> 


■/. 


2R 2 b 

db — U b > 

4 R 2 


and rewrite Equation (A-4) in the form 


M t 

Pa = 77 ?1 0 2 R 2 V <C> — (1 + Z) 6 

P 


or 


P. 


p. o 


c, 0 


<C> (1 + 2)‘ 


where r c 0 = R 2 vj" ! is the collision time for the clouds. 


(A -5) 


50 


4. Conclusion and Outlook 


We now compare these results with the assumptions of the main paper. 

From Equation (A-5) it follows that the annihilation rate is proportional to 
(1 + z) 6 . As we have previously shown, R is almost constant for z > 10. Further- 
more any possible increase of R during the expansion would imply a decrease in 
v. Considering that these quantities cannot change by more than a factor 10 we 
note that their product does not decrease more with decreasing z than the factor 
<cr v > of the main paper. The value of the parameter £ taken as 3 x 10~ 8 in the 
main paper was chosen to fit observational data. The annihilation rate used in 
the main paper is: 

P'a = Vo "p.O <CT V> O + Z > 6 
which can be written in the same form as Equation (A-5) as 


Pi = V< 1 + 2)6 (A ‘ 6) 

a, 0 

where T „ 0 “ ( n p o a v )~ 1 is annihilation time. So that the relation between 
£ and C is obtained from Equations (A-5) and (A- 6) as 

£* = i 


with r „ ~ 10 14 sec and thus 

a , 0 


£ ^ 



(A -7) 


Abell (1961, 1965) and Schatzman (1970) indicate that v is between 10® and 


5.10 8 cm sec -1 and we know that v is decreasing with time. The possible initial 



& 


epenc 


V £~ A 1+ Z 3 a where 0<a<l so that bot^as sump cions lead to 


, .s the^ 
lence of 
Jically more 
;ucn that 

Similar results as 


given in figure 4 to within the uncertainties of the cosmological parameters 
used. 


values of R are given in Table 1 and those which fulfill the relation R < * are 
in the range 10 23 cm to 5 10 23 cm. Evaluation of Equation (A-3) then indicates 
that 10“ 8 < £ < 10“ 6 . The average values of r c 0 are 1.5 10 22 sec and 6.7 10 22 
sec respectively for ?l 0 = 10~ s Mpc" 1 * 3 and 5 10“ 5 Mpc" 3 . Using (A-7) we obtain 
the corresponding values of 

r 3 * 10 -8 for = 10“ 5 Mpc -3 
,lx 10" 8 for 3l 0 = 5 X 10“ 5 Mpc“ 3 


These two values are known with an uncertainty of a factor 5, thus both are 
consistent with the value given in the main paper. We thus conclude that the 
baryon symmetric cosmology of Omnes can lead to a cosmological y-ray 
spectrum of the same form and intensity as that observed. 


*G. Steigman (1969) Has claimed that the y-rays experiments rule out the possibility of equal 
amounts of matter and antimatter in the universe. Wc have already shown that this is not the case. 

More recently Jones and Jones (1970) arrived at the same conclusion as Steigman by consid- 
ering the problem of the collisions. Their main assumption is that half of the mass is annihilated 
(after forgetting the existence of electrons in the annihilation products, and neglecting without 
justification the magnetic effects). It follows from Equations (A- 2 ) and (A-3) that the average 
ratio of the mass annihilated per collision to total mass can be estimated by 


< £ > ~ 


1 m 2 v 2 R 2 

M t 2 x 10 -24 


The order of magnitude of < £ > is ~ 10 -7 . 

Sunyaev and Zeldovich (1970) have considered the question of existence of antimatter on 
large scale in the universe and have shown that, for Q = 1, less than 10 6 of the matter must 
have been annihilated during the period 0 < z < 100. This condition is fulfilled in the model 
discussed here. 

A model considered previously by one of us (Puget 1971) where diffu- 
sion plays the main role in cosmological, annihilation »The results pre- 
sented here indicate that cloud collisions dominate over diffusion in the 
redshift range of significance for -y-ray production so that we now con- 
sider the collision model more physically plausible than the diffusion 
model. 
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Table 1: Neutralization of the Cosmic Plasma 


T 

Peebles (1968) 

This Work 


X 

C (T.x) 

X 

— 

C (T.x) 

4500 

0.92 

1.8 10" 4 

0.96 

3.6 10* 4 

4000 

0.40 

5.9 10 -4 

0.84 

2.2 10~ 3 

3500 

0.072 

2.7 10' 3 

0.64 

7 10‘ 3 

3000 

9.8 10' 3 

2 10* 2 

0.38 

3.2 10~ 2 

2500 

9.2 10" 4 

0.25 

0.14 

0.3 

2000 

1.23 10" 4 

0.96 

0.068 

1 

1500 

5.3 10' 5 

1 

0.06 

1 

1000 



0.049 



Table 2: Size of the Clouds When They Separate 


mm 

n (cm* 3 ) 

Jl (Mpc -3 ) (?I 0 = 5 10" 5 ) 

R (cm) 

\ (cm) 

M (g) 

1000 

3 10 +3 

5 10 4 

5.1 10 22 

3 10 22 

2.8 10 48 

600 

7 10 2 

1.08 10 4 

8.5 10 22 

1.4 10 23 

2.8 10 48 

300 

10 2 

1.35 10 3 

1.7 10 23 

10 24 

2.8 10 48 

150 

10 

1.7 10 2 

3.3 10 23 

10 25 

2.8 10 48 

(Jl 0 = 10“ 5 ) 

1000 

3 10 3 

10 4 

8.6 10 22 

3 10 22 

1.4 10 49 

600 

7 10 2 

2.16 10 3 

1.5 10 23 

1.4 10 23 

1.4 10 49 

300 

10 2 

2.7 10 2 

2.8 10 23 

10 24 

1.4 10 49 

150 

10 

o 

«-( 

CO 

5.6 10 23 

10 25 

1.4 10 49 
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FIGURE CAPTIONS 


Figure 1. Geometry of Intercloud Collision: The interface between the two 
clouds is the y, z plane. The relative motion of the clouds is along the 
x axis. The B field points downward in the z direction and the effective 
current J is in the y direction. The electrons emitted in the annihilation 
and those accelerated by Compton collisions with y -rays loose their energy 
by collision with the photons of the black body radiation generating X-rays. 

Figure 2. The energy-density of the magnetic field at a time t after the be- 
ginning of the collision in either of the two clouds, as a function of x. 
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Figure 1. Cloud-Anticloud Collision, Generation of the Magnetic Field 
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